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SECTION  1 


SEAFARER  CONDUCTIVITY  SURVEY  PROGRAM  SUMMARY 


1 . 1 INTRODUCTION 

A conductivity  survey  has  been  completed  in  the  Upper 
Peninsula  (UP)  region  of  Michigan  with  confirming  measurements 
at  the  Wisconsin  Test  Facility  (WTF) . The  results  define  the 
geological  distributions  within  the  survey  area  and  identify 
average  conductivity  values  for  each  rock  type.  Measurements 
at  the  WTF  provided  a data  base  for  comparative  evaluation  be- 
tween the  scalar  Audio  Magnetotelluric  (AMT)  method  used  in 
this  survey  and  other  methods  such  as  Tensor  AMT,  H/I,  E/H,  and 
Schlumberger . In  addition,  the  WTF  measurement  program  estab- 
lishes a basis  for  scaling  system  current  moment  (IL)  require- 
ments from  an  area  of  known  propagation  performance  to  other 
Michigan  locations. 

The  survey  instrumentation,  support  personnel,  and  com- 
puter analysis  were  provided  by  University  of  Toronto.  Geo- 
logic interpretations  in  this  report  were  made  by  Dr.  D. 
Strangway,  Chairman,  Department  of  Geology,  University  of 
Toronto . 

1 . 2 AUDIO  MAGNETOTELLURIC  (AMT)  TECHNIQUE 

The  AMT  method  measures  the  ratio  of  horizontal  electric 
field  at  the  earth's  surface  to  the  orthogonal  magnetic  field 
to  yield  earth  resistivities  as  a function  of  frequency.  The 
instrumentation  consists  of  a horizontal  dipole  E-field  sensor, 
a coil  H-field  sensor,  a prewhitening  filter  box  and  a vari- 
able frequency  analyzer  receiver.  The  receiver  integrates  the 
E-  and  H-field  signals  and  forms  the  ratio,  from  which  subsur- 
face conductivity  is  inferred.  The  signal  source  is  terrestrial 
spherics  (lightning  discharges).  At  each  measurement  location, 
two  200-foot  wire  antennas  are  laid  out  perpendicular  to  one 
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another,  the  termination  being  a stake  inserted  no  more  than 
12  inches  below  the  surface.  Then  orthogonal  pairs  of  measure- 
ments are  made  at  each  of  ten  predetermined  frequencies  over 
a range  of  three  decades. 

The  output  of  the  receiver  plus  a calibration  factor  is 
equal  to  ten  times  the  logarithm  of  the  apparent  resistivity. 

1 . 3 CONDUCT  OF  SURVEY 

In  the  Michigan  UP  region,  the  low  conductivity  geologi- 
cal formations  tend  to  run  in  a north-south  direction.  For 
this  reason,  the  area  was  sectionalized  into  a grid  configu- 
ration where  measurements  were  made  along  lines  of  latitude 
spaced  approximately  six  miles  apart.  Site  locations  were 
then  evaluated  based  on  a review  of  the  existing  conductivity 
data  base,  geology,  ownership,  and  accessibility.  Sites 
located  in  National  Forests  and  near  population  centers  were 
eliminated.  Selection  criteria  were  based  on  access  and 
ownership  data,  with  State-owned  property  given  highest 
priority.  Figure  1 shows  the  measurement  locations  of  the 
106  sites  in  the  survey  region,  and  Figure  2 illustrates  the 
geology . 

In  Wisconsin,  AMT  conductivity  measurements  were  taken 
along  the  WTF  E/W  antenna  line  at  20  sites  identified  by  the 
existing  antenna  pole  numbers.  The  area  weighted-average 
method  to  define  bulk  resistivities  was  not  suitable  here, 
so  the  resistivities  were  interpolated  between  measurement 
stations  and  then  averaged  at  115  equally  spaced  locations 
along  the  13.1  mile  antenna  line.  Long-wire  AMT  measurements 
were  also  made  utilizing  the  WTF  overhead  antenna  as  the 
E-field  sensor. 

1 . 4 DATA  ANALYSIS 

The  data  are  taken  via  scalar  AMT  measurements  at  sites 
selected  for  their  geological  characteristics  and  proximity  to 
predetermined  lines  of  latitude.  Each  measurement  site  is  as- 
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Conductivity  Measurement  Locations 
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signed  to  a given  line  traverse  to  accommodate  a lateral  pro- 
filing analysis.  The  methodology  of  analysis  is  illustrated 
by  Figure  3,  with  the  major  steps  described  below: 

• Depth  soundinq  curves  are  developed  from  smoothed  data 
which  show  the  actual  field  observations  plotted  sta- 
tion-by-station on  a log  , versus  log  f plot.  This 
method  shows  the  sense  of  resistivity  with  depth  from 
wh  ich  layered  earth  structures  can  be  interpreted. 

• Pseudosections  are  developed  from  a measurement  tra- 
verse of  the  area  and  show  the  resistivity  variations 
as  a function  of  frequency  and  distance.  The  develop- 
ment of  pseudosections  is  an  important  step  in  the 
analysis  in  that  it  allows  locally  wide  variations  to 
be  subordinated  to  achieve  a resistivity  profile  and 
countour  maps  that  can  be  reasonably  integrated  for 
ELF  antenna  performance  predictions. 

From  this  basic  analysis,  the  survey  area  can  be  character- 
ized by  the  following  data  interpretations  and  illustrations: 

a . Lines  Analysis 

Lateral  profiles  are  developed  which  show  resistivity 
gradients  along  the  measurement  traverse  (residuals) , 
layered  earth  interpretations  and  anisotropy  indications 

b.  Geology  Analysis 

The  geologic  distributions  within  the  survey  area  can  be 
used  as  a basis  to  represent  the  total  area  if  the  geolo 
gy  is  well  charted  and  the  resistivities  of  the  various 
rock  types  are  known.  This  requires  many  measurements 
to  be  taken  in  each  rock  type. 

c . Conductivity  Countour  Map 

Contours  of  constant  conductivity  for  the  survey  area 
are  generated  from  the  data  along  a given  measurement 
traverse  and  interpolation  between  measurement  locations 
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along  each  traverse  and  between  lines  of  latitude. 

In  areas  where  data  are  sparse,  geology  maps  are  used 
as  an  aid  in  contouring. 

The  average  resistivities  for  the  survey  area  can  be  com- 
puted for  each  of  these  three  approaches.  Letting  A-j  be  the 
individual  line  traverse  lengths  or  subordinated  areas  associ- 
ated with  each  method,  and  be  the  respective  resistivity, 
the  area  weighted  average  is 


P area  = 

The  area  weighted  averages  associated  with  line  analysis 
and  geology  are  useful  in  bounding  the  analysis  problem  in 
that  one  technique  gives  a one-dimens iona 1 profile  along  ten 
lines  of  latitude  with  no  definition  between  lines  (6  miles 
apart) , while  the  geology  method  characterizes  the  survey  area 
by  rock  type  based  on  the  heuristic  assumption  that  each  rock 
type  and  its  associated  areas  can  be  correctly  identified  (geo- 
logic noise) . The  accuracy  of  each  method  in  characterising 
the  survey  area  is  an  unknown  quantity  since  the  analysis  is 
really  a two-dimensional  problem  and  the  degree  of  geologic 
"noise"  is  always  subject  to  question.  However,  the  combina- 
tion of  these  approaches,  heavily  weighted  by  the  line  analysis 
pseudosections,  can  give  the  best  available  answer  from  measure- 
ment data  in  achieving  the  main  survey  objective,  i.e.,  a two- 
dimensional  conductivity  contour  map  and  its  associated  weight- 
ed average  conductivities. 

Knowledge  of  the  effective  resistivity  pe  is  necessary  to 
determine  the  ELF  excitation  factor  in  order  to  predict  propa- 
gation effects  from  a given  site  and  thereby  size  the  transmit- 
ter array.  The  excitation  factor  is  given  by  the  following 
equation : 


n 

l Ai 
1 


n 

l 


Ai 


ohm-meters . 
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E = 


1 


/h2o (c/v) 


/pe 

► h ‘ c/v 


Estimation  of  the  values  of  the  excitation  factor  and  atmos- 
pheric attenuation  rate  (i.e.,  E and  a)  for  the  purposes  of 
SEAFARER  sizing  has  been  conducted  largely  with  data  from  one 
controlled  source,  the  WTF.  The  estimates  presently  in  use 
are  the  result  of  an  analysis  conducted  by  NUSC  using  a model 
which  shows  that  E and  the  attenuation  rate,  a,  are  coupled, 
which  implies  that  some  knowledge  of  pe  is  required.  Because 
the  GTE  Sylvania  (1976)  AMT  measurement  at  the  WTF  yielded  a 
conductivity  of  0.26  millimhos/meter  under  the  East-West  an- 
tenna as  opposed  to  the  NUSC  (1973)  H/I  measurement  of  0.32 
mil 1 imho/meter , a correction  factor  is  required.  Consequent- 
ly, for  the  purpose  of  sizing  systems,  the  effective  conduc- 
tivity is  taken  to  be 


o 


e 


3.2 

2.6 


1 

►"’area 


This  approach  decreases  the  excitation  factor  used  for  sizing 
by  0.9  dB  from  that  which  would  be  obtained  by  direct  applica- 
tion of  the  1976  AMT  measurements  and  should  therefore  be  con- 
servative. This  scaling  factor  has  not  been  applied  to  the 
data  presented  in  this  report. 

1. 5 SUMMARY  OF  RESULTS 

The  following  sections  summarize  the  survey  results  for 
Michigan  and  the  WTF. 

1.5.1  Michigan  Upper  Peninsula 

The  rocks  of  the  central  UP  region  have  low  conductivity, 
reflecting  the  presence  of  massive  pre-Cambrian  gneisses.  To 
the  east,  the  conductivity  increases  somewhat,  indicating  the 
presence  of  conductive  Cambrian  and  Ordovician  sediments.  This 
limestone  and  sandstone  cover  is  relatively  thin,  since  the  con- 
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ductivities  still  remain  fairly  low  at  2 - 5 x 10“^  mhos/meter. 

At  a depth  of  a few  kilometers  the  conductivity  increases  sharp- 
ly to  values  around  .01  mhos/meters. 

These  patterns  are  illustrated  by  the  95  Hz  contour  map 
shown  in  Figure  4.  The  contour  map  is  based  on  averaged,  smooth- 
ed data  for  106  stations.  Larger  scale  maps  are  to  be  found  in- 
side the  back  cover  of  this  volume.  On  the  95  Hz  map,  a set  of 
observations  taken  by  NU SC  at  75  are  also  shown  and  included  in 
the  contouring.  These  observations  are  in  general  agreement  with 
ours  and  completely  support  the  above  conclusions.  Tables  I,  II, 
and  III  summarize  the  measurement  results  of  this  program. 

1.5.2  Wisconsin  Test  Facility 

AMT  conductivity  measurements  were  taken  along  the  WTF  E/W 
antenna  line  at  20  locations  (200'  E-field  sensors).  Profiles 
of  resistivity  versus  frequency  were  plotted,  and  the  data  were 
smoothed  using  a third-order  polynomial  fit.  Pseudosections  were 
generated  with  the  smoothed  data  and  contoured  by  interpolation 
due  to  uneven  station  spacing.  Data  between  measurement  stations 
were  derived  from  linear  interpolation  of  the  square-root  of  re- 
sistivity. A resistivity  profile  versus  distance  was  generated 
with  data  points  every  600  feet  along  a 13.1  mile  E/W  antenna 
line.  The  115  data  points  were  equally  weighted  and  averaged  to 
determine  the  effective  conductivity  as  seen  by  each  half  and 
full  E/W  antenna  line. 

Half  and  full  antenna  line  AMT  measurements  were  also  made 
utilizing  the  overhead  WTF  lines  as  E-field  sensors.  Over  250 
measurements  were  taken  during  the  period  of  Octover  1975  - 
January  1976  utilizing  morning  and  afternoon  spherics.  The  WTF 
antenna  configuration  is  schematically  illustrated  by  Figure  5. 

The  E/W  antenna  long  and  short  line  test  results  indicate 
that  the  geological  layered-earth  structure  consists  of  three  to 
four  layers;  low  resistivity  at  the  lower  frequencies  is  typical- 
ly about  1000  ohm-meters.  This  is  overlain  by  a high  resistivity 
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Figure  4.  Conductivity  Contours  - 95  Hz. 


TABLE  1 


MICHIGAN  CONDUCTIVITY  SUMMARY 
BASED  ON  E/W  LATITUDE  PROFILES 
(Frequency  - 95  Hz) 


LINE 

EXTENT  OF 
INTERPOLATION 
(MILES) 

P 

OHM-METERS 

a 

MHOS/METERS  x 10-4 

A 

44 

31 

323.00 

B 

56 

355 

28.00 

C 

68 

2884 

3.47 

D 

46 

2691 

3.72 

E 

95 

3802 

2.63 

F 

94 

2884 

3.47 

G 

98 

2630 

3.80 

H 

94 

1660 

6.02 

I 

86 

955 

10.50 

J 

71 

646 

15.50 

AREA  WEIGHTED  AVERAGE 

- 5.67  x 10-4  MHOS/METER 

TABLE  II 


MICHIGAN  CONDUCTIVITY  SUMMARY 
BASED  ON  BEDROCK  GEOLOGY  MAP 
(Frequency  - 95  Hz) 


AREA 
SQ  MILES 

0 

OHM-METERS 

a 

MHOS/METER  x 10~4 

LOWER  6.  MIDDLE  PRECAMBRIAN 

GN,  GR 

1244 

5747 

1.74 

GKM,  GKA , GB,  GH , GKH , H, 
DR 

900 

1020 

9.8 

GK2,  G2 

468 

893 

11.2 

S2  , SM 

324 

238 

42.0 

PALEOZOIC 

OAT,  ELS,  EM,  EJ 

456 

1538 

6.5 

AREA  WEIGHTED  AVERAGE  - 4.48  x 10-4  MHOS/METER 

12 


TABLE  III 


MICHIGAN  CONDUCTIVITY  SUMMARY 
BASED  ON  UNIVERSITY  OF  TORONTO  LOG  p CONTOUR  MAP 
(Frequency  - 95  Hz) 


CONTOUR 
RANGE 
o x 10-4 

AREA 
SQ  MILES 

P 

OHM-METERS 

o 

MHOS/METERS  x 10"4 

a < 2.1 

599 

10,129 

0.987 

2.1  - 10 

1479 

2,197 

4 .55 

10  - 100 

920 

435 

23.0 

a > 100 

431 

70 

143.0 

AREA  WEIGHTED  AVERAGE*  - 5.06  x 

10~4  MHOS/METER 

1 

♦Conductivity  used  for  system  sizing. 
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* x x x x OVERHEAD  TELEPHONE  CABLE 
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Figure  5.  WTF  Antenna  Configuration 
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layer  with  peak  values  in  the  range  of  7,000  to  10,000  ohm- 
meters.  Nearer  the  surface,  the  resistivity  appears  to  become 
lower  with  a hint  that  locally  there  is  a high  resistivity 
surface  layer. 

This  leads  to  the  following  one-dimensional  layered  earth 
structure  interpretation  for  the  WTF: 

Layer 

1 Surface  Strongly  anisotropic  - a few  meters 

thick 

2 P2  % 1700  ohm-meter;  h2  % 600  meters 

3  P3  % 6600  ohm-meter;  113  % 4600  meters 

4 p.  « 1000  ohm-meter 

Recent  (1976)  long  line  AMT  measurements  along  the  N/S  an- 
tenna add  a second  dimension  to  the  above  layered  earth  inter- 
pretation. The  results  show  that  areas  to  the  north  and  west 
exhibit  low  conductivities.  The  GTE  Sylvania  Schlumberger  (1972) 
data,  combined  with  the  E/W  residual  plots  and  measurements, 
suggest  a general  E/W  geologic  fracturing  in  the  bedrock.  A 
simplistic  model  would  consist  of  a shear  line  intersecting  the 
north  leg  1500  meters  from  the  center  and  extending  to  the  west 
leg  5600  meters  from  the  center.  Everything  north  of  the  shear 
line  would  be  highly  conductive.  Anisotropy  near  the  center  of 
the  array  could  be  as  high  as  2. 

The  WTF  survey  emphasis  was  on  measurements  along  the  E/W 
antenna  line.  Data  were  taken  at  20  locations  to  laterally  pro- 
file the  area,  and  full  antenna  measurements  were  made  to  deter- 
mine the  overall  antenna  effective  conductivity.  The  results 
of  these  measurements  and  analysis  indicate  that  the  effective 
conductivity  along  the  E/W  antenna  is  2.6  x 10-^  mhos/meter  at 
80  Hz.  The  results  of  this  survey  and  various  other  WTF  studies 
are  summarized  in  Table  IV. 


SECTION  2 


AUDIO-FREQUENCY  MAGNETOTEELUR1C  (AMT)  THEORY 


2 . 1 INTRODUCTION 

The  magnetotelluric  method  for  determining  subsurface 
electrical  conductivity  structure  was  first  recognized  in 
1950.  Since  then,  most  magnetotellur ic  studies  have 
utilized  low  frequencies  to  investigate  ueep-seated 
features  such  as  sedimentary  basins  or  the  earth's  crust 
and  mantle.  However,  Cagniard  (1953)  explicitly  outlined 
the  extension  of  this  technique  to  shallow  exploration  Ly 
using  higher  frequencies. 

During  the  early  summer  of  1963,  initial  field 
observations  were  made  by  Kennecott  Copper  Corporation  to 
determine  the  feasibility  of  determining  ground  resistivity 
by  magnetotel lur ic  measurements  in  the  audio-frequency  range. 
The  field  program  was  carried  out  in  conjunction  with  DECO 
Electronics,  Boulder,  Colorado,  now  a subsidiary  of 
Westinghouse . Sufficient  encouragement  was  realized  then  to 
justify  Kennecott' s development  of  equipment  specifically 
for  the  purpose. 

For  mineral  exploration  where  resistivity  data  alone 
are  diagnostic,  AMT  has  been  found  to  have  both  advantages 
and  disadvantages  compared  with  conventional  electrical 
methods . 

2 . 2 PRINCIPLE  OF  TEE  METHOD 

The  magnetotelluric  method  has  been  described  by 
Cagniard  (1953),  Wait  (1962),  Vozoff  (1972),  and  others. 
Basically,  the  electromagnetic  impedance  - the  ratio  of  the 
horizontal  electric  field  (E)  in  the  ground  to  the  orthogonal 
horizontal  magnetic  field  (H)  - is  measured  at  a number  of 


frequencies  to  yield  earth  resistivities  as  a function  of 
frequency,  resulting  in  a form  of  depth  sounding.  Operation- 
ally the  technique  is  quite  simple  in  that  depth  sounding  is 
conducted  by  making  measurements  at  only  one  location. 

The  magnetotellur ic  field  can  be  described  in  terms  of 
either  a propagating  electromagnetic  wave  interacting  with 
the  earth's  surface,  or  an  inductive  phenomenon  in  which 
telluric  currents  are  induced  by  the  fluctuating  geomagnetic 
field.  Whereas  the  incident  horizontal  magnetic  field  is 
roughly  doubled  at  the  surface  and  uniform  in  the  absence  of 
lateral  contrasts  in  resistivity,  the  electric  field  is 
directly  dependent  upon  the  earth's  conductivity  structure. 

The  depth  of  sounding  can  be  related  to  frequency  by 
using  the  concept  of  the  skin  depth  given  as 
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where 

6 is  skin  depth  in  meters 
f is  frequency  in  hertz 
o is  conductivity  in  mhos/m 
V is  magnetic  permeability  in  henries/m 

Plots  of  skin  depths  as  a function  of  resistivity  and 
frequency  are  shown  in  Figure  6.  This  figure  shows  that  the 
audio-f requency  range  of  (10  -Hz  to  10  kHz)  covers  the  desired 
depth  range.  As  such,  the  AMT  method  can  be  considered  as 
a supplement  to  other  methods  for  determining  near-surface 
resistivities  such  as  conventional  DC  or  AC  resistivity 
surveys . 

The  apparent  resistivity  can  be  determined  simply  from 
the  measured  impedance  (E/H)  by 

1.26  X 105  X (E/H) 2 , 
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Figure  6 . Skin  Depth  Curves 
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where  E is  in  volts/m,  and  H in  ampere-turns/m.  Data  can 
then  be  plotted  on  log-log  paper  to  form  the  standard 
magnetotelluric  sounding  curve. 

2 • 3 SOURCE  FIELDS 

The  problem  of  source  fields  in  the  audio-frequency 
range  has  been  considered  by  many  workers  in  studies  of 
r^.dio  propagation.  Interest  in  natural  noise  fields  is 
therefore  considerable,  and  several  extensive  studies  have 
been  reported.  The  audio- frequency  magnetotelluric  tech- 
nique commonly  used  in  geophysical  exploration  also  makes 
use  of  natural  audio-frequency  signals  at  150  Hz  and  500  Hz. 
It  is  well  known  and  understood  that  the  main  source  of  nat- 
ural noise  is  thunderstorms,  which  are  extensively  concen- 
trated in  the  tropics  and  tend  to  peak  in  their  activity  in 
the  early  afternoon,  local  time.  The  generated  energy  pro- 
pagates around  the  world,  trapped  in  the  waveguide  formed 
between  the  earth's  surface  and  the  ionosphere. 

Full  treatment  of  the  nature  of  this  propagation  has 
been  given  by  various  authors  and  only  a few  pertinent 
details  will  be  given  here.  The  source  itself  contains  a 
wide  spectrum  of  frequencies;  but  due  to  the  characteris t i :s 
of  the  waveguide,  many  of  these  tend  to  be  attenuated.  Thus 
a pulse  at  the  source  spreads  out  geometrically,  as  well  as 
having  certain  preferred  frequencies  or  rejected  frequencies. 
Among  the  better  known  of  these  are  the  Schumann  resonances 
at  8,  14,  20  and  25  Hz  in  which  fairly  strong  energy  peaks 
are  observed.  The  waveguide  also  has  a strong  absorption 
at  about  2 kHz. 

Figure  7 shows  horizontal  magnetic  field  strength 
curves  for  morning  and  afternoon  on  a typical  summer  day  in 
the  western  United  States  and  Canada.  Note  that  in  all  cases 
the  amplitude  decreases  with  increasing  frequency  up  to  about 
2 kHz,  but  that  this  effect  is  much  greater  in  the  morning, 


20 


0 


FREQUENCY  (MU 


□ 

& 

A 


UTAH.  MARCH.  H( 

A R ZONA  JULY,  96S  A*TfRNOON 
AAiJONA.  JULY,  'MVMO^N'tiG 

BRiTiSh  COLUMB  A.  CANADA.  AwGUST  '966.  AY  T ( •NuON 
•BlTiSH  COLUMBIA.  CANADA.  AUGUST  *66  MOftNM* 


V 


Figure  7.  Typical  Magnetic  Field  Strength  Spectra 


due  to  the  fact  that  the  sources  are  much  farther  away  and 
the  waveguide  attenuation  at  2-3  kHz  is  relativ  ly  much 
greater.  In  the  afternoon,  local  thunderstorm  activity  is 
pri  enl  tnd  the  decrease  at  2 kHz  I Less  Important. 

Oth<  • rw Lse,  the  absolute  field  trength  evel  Is  dependent 
on  the  proximity  to  sources  so  that  in  general: 

a.  Levels  are  higher  in  lower  latitudes; 

b.  Levels  are  higher  in  th<  ift<  rnooi  I 
morning;  and 

c.  Levels  are  higher  in  summer  than  in  winter. 

2 . 4 INTERPRETATION 

The  simplest  method  of  interpreting  magnetote 1 1 ur ic 
data  is  the  comparison  of  the  sounding  curves  obtained  in 
the  field  with  theoretical  curves.  Several  workers  have 
published  theoretical  apparent  resistivity  versus  frequency 
curves  for  two-  and  three-layer  cases,  and  the  calculation 
for  multi-layer  cases  is  quite  simple.  Typical  curves  for 
a two-laye,r  configuration  are  given  in  Figure  f!  . Because 
of  the  law  of  electromagnetic  similitude,  these  layered-media 
curves  can  be  used  for  any  combination  of  layer  thickness  and 
resistivity.  Specifically,  if  we  consider  any  configuration, 
and  multiply  all  lengths  by  and  all  resistivities  by  K , 
the  original  curve  is  correct  if  all  frequencies  are  multi- 
plied by  the  factor 


This  operation  essentially  keeps  constant  the  electromagnetic 
response  parameter  kd  (or  kd'  = i„..  d’  , where  d is  a repre- 
sentative length) . 

Theoretical  curves  to  show  the  effect  of  a layer 
sandwiched  between  an  upper  layer  and  a lower  half-space 
are  given  in  Figure  9.  Essentially,  it  is  net  possible  to 
detect  a conducting  middle  layer  which  is  1/100  the  thickness 
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of  the  upper  layer.  On  the  other  hand,  if  the  thickness 
of  the  middle  layer  is  as  little  as  1/30  that  of  the  upper 
layer  and  it  is  a very  good  conductor,  clear  indications 
should  In • present . Layers  oi  r<  , 

are  hard  to  detect;  and  it  is  not  until  th>  thickness  of 
the  middle  layer  is  two  or  three  times  that  of  the  upper 
layer  that  it  becomes  detectable  as  a separate  layer. 
Limitations  in  measurements  thus  lead  to  the  fact  that  true 
and  accurate  depth  sounding  is  very  difficult.  This 
particular  method  is,  however,  as  are  other  inductive 
techniques,  well-suited  to  searching  for  low-resistivity 
zones . 

Fundamentally,  electrical  current  tends  to  flow  into  a 
medium  of  high  conduct  Lvity;  th<  for<  , current  Induced  bj 
the  geomagnetic  field  may  flow  in  a direction  controlled  by 
the  local  geology  rather  tiian  in  a perpendicular  direction, 
as  expected  when  no  lateral  resistivity  contrast  is  present. 
Moreover,  lateral  resistivity  vari  tions  result  in  a vertical 
magnetic  field  component,  as  used  in  AFMAG . Therefore,  field 
complications  arise  because  the  measured  impedance  at  the 
surface  of  the  earth  is  no  longer  independent  of  either  (1) 
the  orientation  of  the  orthogonal  fields  measured,  or  (2)  the 
polarization  of  the  incoming  source  field.  As  a result,  the 
measured  apparent  resistivity  in  the  vicinity  of  lateral 
resistivity  contrasts  can  have  very  poor  repeatability  in 
time,  especially  ver  period  l . irc< 

thunderstorms  may  have  changed  position.  Similarly,  two 
measurements  using  field  recording  orientations  rotated  90 
degrees  can  show  drastic  differences  in  the  observed  apparent 
resistivity  curves. 

Because  the  magnetotel lur ic  impedance  is  a tensor 
quantity  for  a two-dimensional  earth  resistivity  structure, 
simultaneous  phase-sensitive  measurement  of  both  orthogonal 
components  of  the  electric  and  magnetic  fields  is  required 
for  a complete  description.  Broad-band,  four-channel 
recording  and  subsequent  Fourier  analysis  is  commonly 
utilized  in  low-frequency  magnetotel lur ic  studies,  as 


discussed  by  Cantwell  (I960),  Bostick  and  Sm.th  (1962),  and 
Swift  (1967),  among  others. 
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The  light , field-j  rtal  ised  luri 

survey  was  designed  to  measure  on)\  scalar  resistivities; 
thus,  strictly  valid  dut  t art'  obtained  only  when  the  measur- 
ng  orientati  t!  I 1 r<  tivit;  tructure. 

By  making  closely  , , the  i tioi 

lateral  contacts  can  be  exceedingly  well  defined  and  those 
regions  where  anisotropies  exist  can  be  clearly  delineated. 

With  the  AMT  apparatus  aligned  with  the  resistivity 
structure,  it  is  poi  Lbl<  to  ' • isure  the  impedanci 

for  "E-parallel"  >j  " licular"  polarizations,  ten 

which  refer  to  the  electric-source,  field-oriented  parallel 
to  or  across  the  structure,  respectively.  Apparent 
resistivities  obtained  from  these  impedances  can  then  be 
compared  with  theoretical  two-dimensional  curves 
corresponding  to  the  appropriate  source  polarization. 
Theoretical,  apparent  resistivities  for  arbitrary  two-dimen- 
sional resistivity  structures  are  now  available  by  use  of  a 
computer  modeling  technique  developed  by  T.  R.  Madden  of 
M.I.T.  This  procedure  recognizes  the  analogy  between 
Maxwell's  equations  for  a two-dimensional  geometry  and  the 
transmission  line,  and  constructs  an  electrical  network  in 

which  the  circuit  elements  depend  upon  the  resistivity  struc- 
ture and  the  frequency;  and  the  voltage  and  current  relate  to 
the  electric  and  magnetic  field  values.  The  technique  is  dis- 
cussed by  Madden  and  Swift  (1969),  and  is  elaborated  upon  by 
Swift  (1971). 
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SECTION  3 


INSTRUMENTATION 


3.1  E-  (,  H-FIELU  SK:.S"i-:..  AND  IN]  RP1  ENCE  A : R 

The  instrumental  on  used  in  this  survey  consists  of 
two  sensors  (electric  and  magnetic  field),  a prewhitening 
filter  box  and  a variable  frequency  analyzer. 

3.1.1  Eyy  j«__  H-F ie Id  Sense rs 

The  electric  field  sensor  is  simply  a long  cable 
grounded  at  both  ends.  In  general,  either  500  or  200  feet 
of  dipole  was  used  for  the  sensor.  Although  shorter 
lengths  could  have  been  utilized,  it  was  felt  that  the 
antenna  should  be  as  long  as  possible  in  order  to  average 
out  very  local  effects.  In  Michigan  and  Wisconsin,  all 
data  were  taken  using  200-foot  dipoles  due  to  the  problem 
of  laying  out  longer  wires  in  heavy  brush.  A schematic 
diagram  of  the  detection  equipment  is  shown  as  Figure  10. 

The  magnetic  field  sensor  is  i ferrite-cored  coil 
developed  at  the  University  of  Toronto.  It  has  a frequency 
response  that  is  linear  from  10  Hz  t o 10  kHz,  and  contains 
a series  of  separate  segments,  each  with  its  own  preamplif  i.t 
to  reduce  the  total  capacitance  and  thereby  pioi  id<^  s de- 
bandwidth  in  a single  coil. 

Preconditioning  provides  rejection  of  60  Hz  and  pre- 
whitens the  signal  before  presentation  to  the  analyzer. 

This  preconditioning  also  consists  of  band  limiting  the 
signals,  since  large  signals  from  VLF  stations  and  broadcast 
stations  can  saturate  the  AMT  system  while  it  tries  to 
detect  very  small  signals  at  a frequency  of  interest. 

3.1.2  Interference  Analyzer 

The  analyzer  used  was  assembled  by  Kennecott  Copper 
Corporation  and  is  on  long-term  loan  to  the  University  of 
Toronto.  The  unit  contains  a central  oscillator  and  a 
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Figure  10.  Schematic  Diagram  of  Detection  Equipment 


heterodyning  system  to  measure  the  ratio  of  two  signals 
fed  to  it.  The  filtering  characteristics  of  this  system 
are  very  sharp  and  provide  narrow-band  detection  at  the 
center  frequency.  In  addition,  the  skirts  are  very  sharp, 
since  the  system  is  required  to  detect  very  small  signals 
in  the  presence  of  large,  nearby  ones.  Skirt  rejections 
of  10^  or  more  are  inherent.  The  results  are  presented 
in  the  form  of  a logarithmic  ratio  which,  together  with 
the  gain  settings  and  a calibration  factor,  can  be  used  to 
derive  directly  the  value  of  log  resistivity. 

3 . 2 CALIBRATION 

The  AMT  system  is  calibrated  in  a number  of  ways. 

The  coil  itself  is  calibrated  by  laying  a large  single 
loop  antenna  on  the  ground  with  the  coil  vertical  in  the 
middle.  The  single  loop  is  then  energized  at  each  frequency 
to  be  used  with  a known  magnetic  moment,  and  the  coil  output 
measured.  The  rest  of  the  system  is  then  simply  calibrated 
by  measuring  gain  of  the  separate  stages.  As  the  system  mea 
sures  the  ratio  of  E/H  and  the  channels  are  essentially  the 
same,  there  is  no  temperature  effect  to  a first  approximation 
The  coil  is  not  common  to  both  E-  and  H-field  systems  and 
hence  its  temperature  sensitivity  is  important.  However, 
as  the  coil  is  not  particularly  sensitive  to  temperature, 
results  are  nearly  temperature  independent. 

Other  calibration  procedures  used  have  involved 
energizing  the  coil  and  the  E-field  input  simultaneously, 
and  then  calibrating  directly  in  terms  of  the  ratio  measured. 
This  procedure  gives  essentially  the  same  results. 

The  final  check  on  calibration  has  been  to  make  runs  in 
areas  of  known  resistivity  structure.  For  this  purpose,  one 
area  in  Southern  Ontario  and  one  on  the  Salt  Flats  in  Utah 
have  been  used,  the  latter  having  a uniformly  very  low 
resistivity  medium  and,  hence,  providing  an  excellent  cali- 
bration site.  Data  from  this  site  are  shown  in  Section  5 of 
this  volume. 


29 


SECTION  4 


ERROR  ANALYSIS  AND  SYSTEM  SCALING 


4 . I ERROR  ANALYSIS 

An  estimation  of  errors  associated  with  measuring  resistivity 
can  be  analyzed  by  identifying  the  major  sources  of  dispersion 
in  the  result.  One  is  the  normal  dispersion  to  instrumental 
errors  that  can  be  measured  in  the  laboratory,  and  another  is 
dispersion  in  the  resistivity  measurement  itself,  i.e.,  measure- 
ment repeatability. 

The  instrument  can  be  separated  into  the  receiving  coil, 
the  preamplifiers,  and  the  interference  analyzer.  The  contribu- 
tion of  each  to  the  total  error  can  be  measured  and  converted  to 
an  equivalent  resistivity  error.  The  receiving  coil  must  be 
calibrated  with  a known  magnetic  field  at  each  of  the  measurement 
frequencies.  The  accuracy  of  the  calibration  results  in  the 
errors  in  the  resistivity  of  2 to  3 percent.  The  variation  in 
the  true  sensitivity  of  the  coil  deduced  from  calibrations  that 
were  repeated  several  months  later  gives  resistivity  errors  of 
10  to  15  percent. 

The  gain  of  the  preamplifiers  was  measured  before  and 
after  the  Michigan  and  Wisconsin  field  trips.  The  maximum  gain 
errors  were  1.1  dB  resulting  in  resistivity  errors  of  14  percent. 
The  sensitivity  of  the  gain  to  temperature  variations  has  not 
been  measured. 

The  Interference  Analyzer  also  has  measurement  errors  in 
the  value  of  the  ratio  of  voltages  measured  and  in  the  center 
frequency  "F0"  of  its  filter.  The  repeatability  of  the  ratio 
is  within  ±.3  dB,  which  is  equivalent  to  a resistivity  error  of 
17  percent  of  F0.  Center  frequency  F0  is  repeatable  to  within 
+ 2 Hz  (11%  fQ) . 

The  instrumental  and  calibration  errors  give  total 
uncertainty  of  approximately  129  percent  in  resistivity 
readings . 
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The  second  source  of  dispersion  can  be  estimated  by  making 
repeat  soundings  at  one  station  where  the  resistivity  structure 
is  well  known.  Repeatability  measurements  were  made  on  the 
Bonneville  Salt  Flats  in  Utah  where  the  resistivity  is  known 
to  be  in  the  range  .25  to  .32  ohm-meters  essentially 
independent  of  frequency  for  these  particular  measurement 
frequencies . 

The  data  from  the  repeat  readings  at  each  measurement 
frequency  for  all  orientations  of  the  electric  field  sensor 
were  averaged  and  a standard  deviation  calculated,  'in  one 
case, the  square  roots  of  the  resistivities  were  averaged, 
which  is  equivalent  to  averaging  the  electric  fields;  in  the 
other  case, the  resistivities  were  averaged  directly.  The  results 
of  these  calculations  are  shown  in  Figures  11  and  12.  The 
1 kHz  and  2 kHz  data  have  associated  with  them  substantially 
larger  standard  deviations  than  are  observed  at  the  other 
frequencies.  The  average  values  appear  to  come  reasonably 
close  to  other  published  data,  excepting  the  frequencies  just 
mentioned  that  have  anomalously  high  values. 

These  data  provide  a test  of  the  repeatability  of  the 
measurements  at  one  particular  station  during  a six-hour 
measurement  period.  The  different  values  of  resistivity  that 
are  measured  are  not  the  result  of  instrumental  errors  (that 
are  less  than  ±10%  in  this  case) , although  they  incorporate  some 
of  these  error  sources,  but  are  related  more  to  the  fact  that 
during  each  measurement  period  the  instrument  is,  in  fact,  seeing 
different  signals  from  different  sources  where  orthogonal  com- 
ponents of  the  electric  and  magnetic  fields  are  not  necessarily 
related  by  Cagniards’  formula  for  the  apparent  resistivity.  That 
is,  they  may  be  what  is  referred  to  as  uncorrelated  transients 
from  local  sources  and,  hence,  give  "apparent  resistivities"  that 
are  different  from  the  true  apparent  resistivity. 

The  error  analysis  summary  is  shown  in  Table  V.  If  the 
distribution  of  instrumentation  errors  are  all  biased  in  one  di- 
rection, then  the  worst-case  resistivity  error  would  be  29  percent. 
On  the  other  hand,  if  the  error  bias  is  randomly  distributed,  in- 
strumentation errors  would  normally  be  less  than  10  percent. 
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Figure  11. 


Salt  Lake  Test  Data:  Average  Square  Root  of  P 


TAB LB  V 


ERROR  analyst:;  summary 


INSTRUMENTATION 

PERCENT  ERROR  IN  p 

COIL  CALIBRATION: 

Repeatabi 1 ity 

+ 3.0 

True  Sensitivity 

+ 15.0 

PREAMPLIFIER  GAIN 

+ 4.0 

INTERFERENCE  ANALYZER 

+ 7.0 

TOTAL  +29.0 

MEASUREMENT  REPEATABILITY 

Percent  Standard  Deviation 

in  p at  79  Hz 

1 MEASUREMENT 

+ 42.0 

20  AVERAGED  MEASUREMENT  ( WTF ) 

+ 9.39 

106  AVERAGED  MEASUREMENTS 

i 4.08 

(MICHIGAN) 
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The  technique  itself  produces  errors  that  are  larger  than 
the  instrumentation  errors  for  a few  data  samples.  In  the 
Michigan  site  survey,  the  resistivity  data  from  106  measure- 
ment stations  were  averaged,  which  reduces  the  percent  standard 
deviation  in  resistivity  from  42  percent  to  4.08  percent.  In 
the  WTF  survey,  the  percent  standard  deviation  in  resistivity 
is  9.34. 

4.2  SYSTEM  SIZING 

As  indicated  in  section  1.4,  system  sizing  parameters  for 
SEAFARER  have  been  developed  largely  from  data  obtained  from  a 
single  controlled  source,  the  WTF.  Sizing  of  systems  to  meet 
communication  objectives  for  other  sites  is  basically  a pro- 
cess of  scaling  WTF  performance  as  a function  of  location  and 
•trical  properties  of  the  candidate  site.  Because  the  Nl 
analysis  giving  present  estimates  of  ELF  propagation  para- 
meters used  for  SEAFARER  sizing  recognizes  coupling  between 
excitation  factor  E and  attenuation  rate  i and  shows  that  . is 
a slowly  varying  function  of  E,  it  has  been  deemed  appropriate 
to  employ  the  expression 

3.2  1 


developed  in  section  1.4.  This  scaling  restores  the  self- 
consistency  of  the  propagation  parameters  determined  by  NUSC, 
and  yields  s<  wn»  »wl i ■ * tonserval  I ■ factoi 
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SECTION  5 


DISCUSSION  OF  RESULTS 


The  following  sections  discuss  the  AMT  measurement  re- 
sults in  Michigan  and  the  WTF  with  detailed  graphs  and  com- 
puter analysis  contained  in  Appendixes  A and  B,  respectively. 
In  general,  the  discussions  and  graphs  are  subdivided  as 
fol lows : 

1.  Log  Resistivity  vs.  Log  Frequency  with  Polynomial 
Fits 

• North-South  Orientation 

• East-West  Orientation 

• Averaged  Data 

2.  Pseudosections 

• North-South  Orientation 

• East-West  Orientation 

• Averaged  Data 

3.  Residuals 

• North-South  Orientation 

• East-West  Orientation 

• Averaged  Data 

4.  Anisotropy  (TE/TM  Ratio) 

5.  Measurement  Locations,  Latitude  and  Longitude 

5.1  MICHIGAN  UPPER  PENINSULA  RESULTS 

5.1.1  Depth  Sounding  Curves  - Log  . vs.  Log  f 

The  first  set  of  data  shown  in  Appendix  A.l  is  the 
actual  field  observations  plotted  station  by  station  on  a 
log  resistivity  - log  frequency  plot.  This  system  shows  the 
sense  of  resistivity  with  depth  in  a general  way  in  that  de- 
creasing frequency  means  increasing  depth.  There  are  106 
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st  it  i ' r i s of  late.  Results  were  obtained  using  200-foot 
elect  iodi  arrays  (end  to  end)  oriented  both  N/S  and  E/W.  The 
: ua  1 • n each  plot  ire  the  actual  observations.  In  genera)  , 

it  : ■ • • i • ;i":;  '■  in  ; ■ ■ \ irirn  liately.  Th< 

resist  • wit  frt  juency  foi  most  stations.  At  a sig- 

nificant number  of  stations,  there  is  also  a low  at  approxi- 
mately 2-4  kHz.  This  low  may  or  may  not  represent  true  ground 
resistivity  since  the  signal  strengths  become  very  weak  at 

There  is  also  a gen- 
eral smoothing  of  the  curves  when  both  E/W  and  N/S  observa- 
tions are  combined  into  one. 

The  curves  shown  on  these  plots  represent  a smootning 
of  the  observations  using  a third  order  polynomial.  Such 
a polynomial  roughly  smooths  data  to  the  level  of  resolu- 
tion that  can  be  expected  for  magnetote 1 lur i c observations. 
Subsequent  plots,  therefore,  use  the  data  derived  by  smooth- 
ing as  their  basis. 

Stations  marked  with  a K represent  repeat  observations 
or  poor  data  and  are  simply  included  here  for  the  record  i :d 
not  used  for  analysis. 

5.1.2  Pseudo sect  ions 

The  pseudosections  shown  in  Appendix  section  A. 2 
portray  the  resistivity  distribution  in  a two-dimensional 
array.  Lower  frequenev  represents  greater  depths. 

The  data  are  contoured  using  interpolation  between  stations 
at  two-minute  intervals.  This  interpolation  is  based  on  the 
square  root  of  the  resistivity  and,  therefore,  simulates  the 
actual  observations  based  on  the  electrical  field.  As  one 
moves  from  line  to  line,  some  consistent  pattern  can  be  seen. 
First,  line  A is  highly  conductive  compared  to  the  other  lines. 
As  we  will  see,  this  is  due  to  the  fact  that  most  of  Lire  A is 
on  Precambrian  sediments  or  on  slates  which  are  relatively 
conductive.  Further  south  the  resistivities  are  generally 
higher  in  the  main  Precambrian  gneiss  and  granite  areas. 
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A number  of  sharp  lateral  contrasts  exists,  but  the 
decrease  with  depth  shows  up  consistently  on  all  pseudo- 
sections . 

5 . 1.3  Residua  1 Plots 

^he  residual  plots  are  a form  of  data  presentation  de- 
signed to  emphasize  the  differences  along  lines.  The  pro- 
files are  based  on  interpolated  data  and  are  shown  in 
Appendix  A. 3.  On  the  basis  of  these  interpolations,  it  is 
possible  to  derive  a mean  resistivity  for  each  line  at  each 
frequency.  This  value  is  based  on  th<  Length  thi  Lin 
which  observations  were  made  and  was  determined  by  averaging 
» . at  each  two-minute  interval.  These  values  are  tabulated 
as  a function  of  frequency  for  each  line  traverse  in  Table  VI. 
Table  VII  summarizes  the  area  weighted  averages  for  all  lines 
as  a function  of  frequency. 

In  the  depth  sounding  and  residual  plots,  the  pattern 
of  lower  resistivities  at  low  frequencies  and  the  sense  of 
increasing  resistivity  in  the  middle  lines  are  seen.  One  -an 
also  quite  clearly  see  the  presence  of  an  additional  low  re- 
sistivity layer  at  the  hiqher  frequencies  in  the  range  of 
1 to  4 kHz. 

These  values  have  in  turn  been  used  to  derive  an  over  ill 
weighted  curve  for  the  region.  This  curve  is  weiqhted 
in  terms  of  the  length  of  the  lines  and  hence  represents 
an  overall  weiqhted  average  for  the  region  surveyed,  the 
weighting  being  in  terms  of  the  area  sampled.  The  line 
average  and  area  weighted  resistivity  plots  are  shown  in 
Figures  13,  14,  and  15. 

A preliminary  attempt  has  been  made  to  interpret  this 
as  a two-layer  structure.  A rough  approximation  is  for  a top 
layer  two  kilometers  thick  with  a resistivity  of  about  2000 
ohm-meters.  This  overlies  a very  good  conductor  which  has  a 
resistivity  less  than  50  ohm-meters.  While  there  are 
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significant  variations  from  place  to  place,  this  overall 
weighting  probably  represents  typical,  shallow  Precambri  m 
crust . 

The  residual  profiles  show  the  deviation  from  the  line- 
by-line  means  at  each  frequency.  On  ^he  basis  of  these 
plots,  one  can  trace  the  deviations  from  the  average  quite 
well.  i • xant]  Le , the  resistivit  • it  lepth  on  the  west  end 
of  Line  A is  below  the  average. 

On  line  B the  resistivity  at  the  east  end  is  high,  pre- 
cisely where  the  stations  are  located  on  granite  gneiss.  The 
sparse  sampling  on  the  west  end  is  on  slates  or  sediments  and 
shows  low  resistivities. 

Again,  on  Line  C the  granite  gneiss  is  resistive,  the 
slates  conductive.  Stations  S33  and  S34  are  also  relatively 
conductive,  showing  that  the  graywackes  are  conductive.  Station  N12 
is  highly  resistive  and  on  granite.  This  same  pattern  can 
be  seen  on  Line  D;  in  fact,  over  the  whole  area  we  see  that 
the  highs  on  the  residual  plots  essentially  correspond  to  the 
granites  and  the  granite  gneisses  of  the  region. 

Line  E shows  essentially  the  same  pattern  although  sta- 
tion S48  does  not  follow  the  pattern.  At  the  east  end  of  the 
line,  stations  S3,  SI,  S4,  S4B  and  S2  are  all  on  the  Paleo- 
zoic cover.  The  residual  plots  show  the  effect  of  this  quite 
clearly,  the  more  conductive  sediments  thickening  to  the  east 
so  that  at  station  S2  the  data  are  hardly  affected  by  the 
underlying  Prccambrian. 

The  pattern  for  Line  F continues  to  show  the  same  infor- 
mation with  high  resistivities  for  the  major  areas  of  granite 
gneiss.  The  Paleozoic  sedimentary  cover  at  the  east  end  is 
quite  thin  as  seen  both  by  a window  of  Prccambrian  and  by  the 
lack  of  low  resistivities,  suggesting  that  the  sediments  here 
are  quite  thin. 
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Line  G continues  to  have  the  same  pattern  of  highs  over 
the  granite  gneisses  and  lows  over  the  graywackes.  The  Paleo- 
zoic sediments  are  likely  quite  thin  as  they  scarcely  influence 
the  observed  resistivity  values  out  to  station  S9 . 

On  line  H,  the  high  resistivity  regions  shown  again  very 
closely  reflect  the  granite  gneiss  areas.  The  sedimentary 
cover  at  the  east  appears  to  be  substantial  at  stations  S54 
and  S-55  but  not  closer  to  the  Precambrian  contact. 

The  same  situation  is  repeated  on  both  Line  I and  Line  J. 

5.1.4  Anisotropy  TE/TM 

The  final  sets  of  plots  shown  in  Appendix  section  A. 4 
are  the  ratios  of  the  N/S  to  L/W  components  (differences  in 
log  domain)  which  can  be  used  to  study  lateral  variations. 

Where  the  values  differ  substantially  from  each  other,  these 
plots  emphasize  the  differences.  In  general,  the  region  is 
not  highly  anisotropic,  nor  does  it  seem  to  have  an  overall 
regional  effect.  A number  of  the  sharp  changes  do  seem  to 
correspond  roughly  to  sharp  lateral  changes  in  resistivity. 

This  is  not  always  the  case;  but  where  it  is,  it  is  undoubt- 
edly due  to  the  coupling  of  the  source  energy  to  the  later  il 
contact  giving  a large  apparent  anistropy. 

5.1.5  Summary  of  35  Hz  and  95  Hz  Data 

In  order  to  present  a simple  summary  of  the  regional 
variations  in  resistivity,  the  average  value  of  resistivity 
at  35  Hz  and  95  Hz  has  been  plotted  on  two  separate  countour 
maps  (see  Maps  1 and  2,  respectively).  These  data  represent 
the  smoothed  average  values  for  the  106  stations  and  are 
shown  as  log  ; values  which  have  been  contoured  on  an  inter- 
val of  0.33  log  p.  Figure  4 in  section  1 illustrates  four 
ranges  in  conductivity  over  the  survey  region  at  95  Hz.  On 
the  95  Hz  map,  a set  of  observations  taken  by  NUSC  at  75  Hz, 
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are  also  shown  and  included  in  the  contouring.  These  observa- 
tions are  in  general  agreement  with  ours  and  completely  support 
the  foregoing  conclusions.  In  general,  nearby  stations  give 
similar  results,  thus  corroborating  our  data  by  totally  inde- 
pendent observations. 

The  overall  regional  pattern  can  be  summarized  as  follows: 

a.  The  granites  and  granite  gneisses'  are  high  in 
resistivity . 

b.  The  graywackes  and  slates  are  fairly  low  in 
resistivity. 

c.  The  sedimentary  cover  in  the  south-east  part  of 
the  area  is  relatively  low  in  resistivity,  but 

it  appears  to  be  very  thin  and,  hence,  influences 
the  observations  at  only  a few  places. 

5 . 2 WISCONSIN  TEST  FACILITY  (WTF) 

During  the  period  of  October  1975  - January  1976,  AMT 
conductivity  measurements  were  taken  along  the  WTF  E/W  and 
N/S  antenna  lines.  Short  line  measurements  (200’  E-field 
sensors)  were  taken  at  20  locations  along  the  E/W  antenna 
right-of-way  (ROW) . Half  and  full  antenna  measurements  (7/14  mile 
E-field  sensors)  were  also  taken,  utilizing  the  existing  WTF 
overhead  antenna  lines.  The  following  subsections  discuss 
the  results  of  these  measurements  with  particular  emphasis  or 
the  E/W  line  data,  since  this  will  be  used  as  a basis  for 
scaling  system  current  moment  ( I L ) requirements  from  an  area 
(WTF)  of  known  propagation  performance  to  other  locations. 

5.2.1  Short  Electrode  AMT  Results  Along  E/W  An tenna  ROW 

AMT  conductivity  measurements  were  taken  along  the  WTF 
E/W  antenna  line  at  20  locations  (200’  E-field  sensors). 

Profiles  of  resistivity  versus  frequency  were  plotted,  and 
the  data  were  smoothed  using  a third-order  polynomial  fit. 
Pseudosections  were  generated  with  the  smoothed  data  and 
contoured  by  interpolation  due  to  uneven  station  spacing. 
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Data  between  measurement  stations  were  derived  from  linear 
interpolation  of  the  square-root  of  resistivity.  A resis- 
tivity profile  versus  distance  was  generated  with  data  points 
every  600  feet  along  a 13.1  mile  E/W  antenna  line.  The  115 
data  points  were  equally  weighted  and  averaged  to  determine 
the  effective  conductivity  as  seen  by  each  half  and  full  E/W 
antenna  line. 

The  following  paragraphs  discuss  the  measured  results 
from  the  depth  sounding  data  and  the  respective  pseudosections, 
residual  plots,  and  anisotropy  data. 

5.2. 1.1  Log  Resistivity  vs  Log  Frequency  Results 

Appendix  B.l  contains  the  Log  t vs  Log  f plots  for  each 
sounding  made  at  the  WTF . These  consist  of  observations  made 
at  20  stations  along  the  east-west  transmitter  line.  The 
station  numbers  correspond  to  pole  numbers  (which  are  located 
at  roughly  100  meter  intervals).  The  first  letter  on  these 
plots  refer  to  a half  of  the  transmitter  line:  east  (E)  or 

west  (W) . The  next  letter(s)  refers  to  the  electric  field: 

N/S  (N) , E/W  (W) ; or  to  the  average  (AVG) . The  data  points 
shown  refer  to  the  actual  field  calculation  of  the  apparent 
resistivity.  Many  soundings  at  a specific  frequency  were 
repeated.  In  these  cases,  the  value  which  gave  the  smoothest 
curve  in  the  resistivity-frequency  plot  was  chosen.  The  raw 
data  can  be  verified  by  referring  to  the  original  field  data 
on  file  both  at  the  University  of  Toronto  and  GTE  Sylvania. 

It  is  worth  noting  that,  in  general,  the  values  of  re- 
sistivity tend  to  be  low  at  2160  Hz  and  at  the  adjacent 
frequencies  of  864  and  4310.  These  low  values  may  be  related 
to  very  weak  source  fields  encountered  at  the  time  of  obser- 
vation, and  are  believed  to  be  a result  of  problems  in  measur- 
ing these  very  low  signals. 

The  smooth  curve  through  each  set  of  data  represents  a 
least  squares  third-order  polynomial  fit.  The  choice  of  a 
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third-order  polynomi  il  is  Logical,  albeit  somewhat  arbitrary, 
because  it  permits  the  detection  of  four  layers  of  resistivity, 
which  is  about  the  limit  of  detectability  for  the  three  decades 
of  frequency  used.  The  E/W  and  N/S  data  ire  somewhat  "noisier" 
than  the  averaqed  data.  (The  average  was  determined  by  taking 
the  square  root  of  the  resistivity  and  summing  and  squaring.) 
This  calculation  leads  to  the  smoothed  curves  shown  in  the 
appendix . 

At  this  time,  criteria  for  rejectinq  individual  observa- 
tions have  not  been  developed,  and  the  effect  that  individual 
points  have  on  determining  the  polynomials  reman  to  be  studied 
in  the  future.  In  general,  however,  it  is  noted  that  the 
smoothed  polynomial  fits  the  observed  data  points  quite  well, 
especially  at  the  frequencies  of  interest  to  this  project.  (In 
preparing  the  appended  plots  one  extreme  data  point  at  79  Hz  at 
station  EW  112  was  eliminated.) 

The  data  suggests  that  the  structure  consists  of  three 
to  four  layers.  Low  resistivity  tl  frequencies  is 

typically  about  1000  ohm-meters.  This  is  overlain  by  a high 
resistivity  layer  with  peak  values  in  the  ranqe  of  7000  to 
10,000  ohm-meters.  Nearer  the  surface,  the  resistivity  appears 
to  become  lower,  hinting  that  locally  there  is  a high  resis- 
tivity surface  layer.  It  is  pertinent  to  make  a few  further 
observations : 

a.  The  data  were  collected  along  the  WTF  east-west 
line,  and  there  is  some  concern  that  this  mighty 
have  had  some  effect  on  the  observed  results.  In 
general  the  N/S  electric  field  observations  are 
higher  than  the  E/W  observations,  although  there 
are  exceptions.  This  difference  may  be  related 
either  to  a preferred  geologic  qrain  or  to  the 
cable.  One  station  (97W)  was  measured  100  feet 
off  the  line,  and  the  same  pattern  was  observed. 
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which  suggests  that  the  effect  is  geoloqical.  A 
general  E/W  geologic  fracturing  or  jointing,  with 
some  watei  filling  the  sheers,  would  account  for 
these  observations. 

I . At  the  highest  E requenci< 

•ted.  It  seems  possible  that  this  was  a 
result  of  coupling  in  some  way  to  the  overhead 
antenna  cable.  It  is,  of  course,  also  possible 
that  the  local  overburden  has  filled  small  bed- 
rock valleys  that  have  a preferred  orientation. 

5. 2. 1.2  Pseudosections 

Pseudosections  are  an  attempt  to  present  the  variation 
in  resistivity  both  as  a function  of  frequency  and  as  a 
function  of  position  alonq  a lino.  Since  the  soundings  were 
not  made  at  equal  intervals,  an  int  i olation  procedure  was 
developed.  For  the  WTF,  200  meters  (the  distance  between 
56W  and  58W)  was  chosen  as  the  inter;  olation  spacing.  Poly- 
nomial fits  at  each  sounding  point  were  used  as  the  lasis, 
then  interpolated  between  stations  usin;  (resistivity)^  “ . 
This  calculation  method  was  chosen  because  it  most  closely 
represents  an  attempt  to  consider  how  the  electric  field 
would  bi  i . . rhesi  lata  thei  :onl  u l I 

graphic  view  of  the  way  in  which  the  resistivity  varies  lat- 
erally with  frequency.  Appendix  B.2  shows  the  pseudosections 
for  the  E/W  antenna  line. 

Examining  the  N/S  data  along  the  E/W  intenna  line,  one 
sees  immediately  that  at  79  Hz  the  resistivity  is  generally 
in  the  range  of  2,000  to  20,000  ohm-meters.  Tho  iveraqe 
values  are  tal  ted  Cable  VXI1.  Th i p raq 
is  based  on  taking  the  square  root  of  each  resistivity  value 
in  the  interpolated  data  set  and  including  only  one  data 

Lnt  : • fond  the  ends  ©i  th<  bs  rval  Lon  interval.  In  '.".feet 
this  means  that  resistivity  value  is  properl',  weighted  to  con 
sider  the  variations  from  station  to  station. 
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RESISTIVITY  - WTF  (E/W  ANTENNA  LINE) 
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As  noted  earlier,  the  resistivity  drops  to  low  values 
at  low  frequencies  w’th  a low  resistivity  sequence  at  around 
2 kHz.  We  believe  this  t be  the  result  of  weak  fie]  that 
make  the  observations  difficult.  The  pattern  for  the  L/W  ob- 
servations is  similar,  except  that  the  values  tend  to  be  low- 
er due  to  the  anisotropy  already  mentioned. 

5. 2. 1.3  Residual  Plots 

The  residual  plots  shown  in  Appendix  B.3  are  designed 
to  emphasize  the  deviation  from  the  mean  in  the  apparent  re- 
, tivity . f resistivil  termined 

line  have  been  used  as  baseline  values  and  are  tabulated  in 
Table  VIII.  The  dat  i art  intt  rj  <1  ited  on  200  meter  inter- 
vals (note  ticks  at  actual  observation  points),  and  are 
plotted  so  that  the  vertical  scale  is  as  shown  in  units  of 
10  log  ,■  (i  ohm-m)  . An  examination  of  the  profiles  permits 
one  to  quickly  see  the  regions  of  high  and  low  resistivity 
and  how  these  vary  with  frequency.  In  a general  way,  the 
north-south  and  east-west  components  show  quite  similar 
patterns  - for  example,  low  resistivities  at  low  frequencies 
near  station  113W  and  general  lows  near  stations  12W  and  26W. 
Different  positions  for  the  anomalies  at  the  two  polarizations 
are,  however,  quite  clear  and  distinct  in  several  cases.  This 
pattern  is  expected  across  lateral  contrasts  in  resistivities 
and  confirms  the  existence  of  such  contrasts.  A good  example 
of  this  is  at  station  30E  where  the  east-west  component  is 
high  and  the  north-south  component  is  low. 

Takinq  the  square  root  average  of  these  data  results  in 
the  last  figure  of  this  set  in  which  the  deviations  from  the 
mean  are  much  subdued  (note  reduced  scale) , so  that  the  sharp 
effects  at  lateral  contacts  have  been  smoothed.  These  graphic 
displays  are  a useful  way  of  displaying  the  deviations  from 
the  simple  patterns  otherwise  observed. 

5. 2.  1.4  Anisotropy  TE'/TM 

Separate  plots  designed  to  emphasize  the  anisotropy  be- 
tween the  two  observing  directions  have  been  determined  at 


11 


each  frequency  as  a function  of  distance  alonq  the  traverse 
line  and  are  shown  in  Appendix  B.4.  The  data  have  been  in- 
terpolated to  a 200  meter  spacing  and  the  ratio  of  the  two 
components  shown  directly.  Again,  the  vei t ical  scale  is 
logarithmic  so  that  a simple  perspective  is  obtained.  As 
already  noted,  it  is  • 1 1 1 - striking  that  the  north-south 
resistivity  is  consistently  hiqher  than  the  east-west  resis- 
tivity. This  implies  a strong  anisotropy  to  the  whole  region, 
an  effect  which  can  be  attributed  largely  to  the  presence  of 
east-west  oriented  fracturing,  perhaps  controlled  by  a few 
major  shear  zones  in  that  general  direction,  it  is  possible, 
of  cour  , that  sonu  this  i:n  s>  ' r ' ; i:;  d .i  - ■ i - ) : v iv- 

sence  of  the  transmitter  cable,  which  seems  to  be  unlikely, 
except  at  the  two  highest  frequencies  where  the  discrepancy 
is  considerable.  Only  at  st  ti ion  26W  is  there  a general  local 
reversal  of  the  anisotropy  pattern  of  high  values  in  a north- 
south  direction  and  low  values  in  an  east-west  direction. 

Again,  this  discrepancy  appears  to  be  a result  of  sharp  lat- 
eral changes  as  already  seen  in  the  residual  plots  (see 
Appendix  B. 3 ) . 

.2.2  Hal  f_  and Fu  1 1 Antenna  AMT  Mi  lasu  r»  r 

During  the  period  of  October  1975  - January  1976,  AMT 
measurements  were  taken  at  the  WTF  utilizing  the  E/W  and  N/S 
antenna  lines  as  E-fiold  sensors.  Over  250  measurements  were 
taken  during  this  period  utilizing  morning  and  afternoon 
spherics.  The  E-field  sensor  configuration  is  tabulated  in 
" ibl  v t md  • h<  w'i’i  'erhead  ant  . :hei  I ?a 

illustrated  by  Figure  5. 

The  measurement  data  were  weighted  equally  for  frequencies 
70,  85,  and  95  Hz. 
usual  manner  as  follows: 


data  points. 


0 = 


1 n • 

5 451 


number  of 
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TABLE  IX 


AMT  E-FIELD  SENSOR  CONFIGURATION 


E-FIELD 

SENSOR 

CONFIGURATION 

MEASUREMENT 

LOCATION 

East  Leq 

P112E  to  WTF 

P112E 

West  Leq 

P112W  to  WTF 

P112W 

E/W  Antenna 

P112E  to  West  Terminal  Ground 

P112E 

W/E  Antenna 

P112W  to  East  Terminal  Ground 

P112W 

North  Leq 

P100N  to  WTF 

P100N 

South  Leq 

P74S  to  WTF 

P74S 

N/S  Antenna 

P100N  to  South  Terminal  Ground 

P100N 

S/N  Antenna 


P74S  to  North  Terminal  Ground 

P74S 
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Data  were  not  smoothed,  since  a large  percentage  of  the  measure- 
ments were  taken  over  a limeted  frequency  range  between  35  and 
213  Hz.  The  measurement  results  are  summarized  in  Table  X 
with  he  mean  and  standard  deviations  listed  in  Table  XI. 

The  existing  N/S  overhead  antenna  has  in  its  near  proximity 
several  cables  which  could  be  considered  parasitic  elements. 
Figure  14  illustrates  the  situation  where  a faulted  buried  line 
is  within  25  meters  of  the  overhead  antenna  over  a 7-mile  length 
of  the  north  line.  In  addition,  a recent  rerouting  of  telephone 
cable  from  Hwy  GG  to  the  north  leg  ROW  places  a parasitic 
element  within  10  meters  of  the  overhead  line  over  a length  of 
2-1/2  miles.  Such  proximity  of  these  cables  is  bound  to  in- 
fluence the  E/H  wave  impedance  measurement  as  well  as  any  WTF 
propagation  or  pattern  steering  experiments.  Recent  NUSC 
pattern  steering  measurements  bear  out  the  fact  that  the  WTF 
antenna  has  changed  since  1973.  (Private  Communication, 

P.  Bannister,  NUSC)  Phasing  of  the  various  antenna  elements 
are  different  by  20°  to  30°  from  previously  measured  a. i-  i- 
thal  patterns. 

The  south  overhead  antenna  does  not  seem  to  be  affc'ted 
by  the  faulted  buried  cable  near  the  overhead  line.  Thi  is 
probably  due  to  the  fact  that  the  buried  cable  is  >nly  n ‘h-._ 

ROW  over  5 percent  of  its  length  and  has  an  average  sparine 
of  0.26. 

The  measurement  results  for  the  conductivity  under  the 
south  antenna  leg  is  consistent  with  previous  studies,  while 
ifj  and  ' N/s  values  are  questionable,  suggesting  parasitic 
contamination.  Measurement  errors  associated  with  contamina- 
tion could  be  caused  by  one  of  the  following  factors: 

a.  It  has  been  shown  from  modeling  and  measurement 
verification  that  a grounded  long  wire  parasitic 
line  (Zp  = 100  ohms)  located  within  0.016  (40  meters) 
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Possible-  parasitic  contami  nat  ion. 


tabu:  xi 


STATISTICAL  ANALYSIS  OK  Till  W i F MEASUREMENT  DATA 


ANTENNA 

CONFIGURATION 

X 

MEAN 

OHM- METERS 

STANDARD  DEVIATION 
OHM-METERS 

North  Lcq 

1556 

ui 

U1 

Sout  h Leq 

7802 

2625 

East.  Leq 

4392 

1945 

West  Leg 

3721 

1893 

N/S  Predicted 

4679 

2687 

N/S  Measured 

2497 

1086 

E/W  Predicted 

4056 

2714 

E/W  Measured 

4065 

1305 

E:  ; ; ■ I i ■ ' i -i  iard  ■ Lai  tnu  rroally 

distributed  independent  half  lea  distributions. 


I 

I 

I 

I 
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of  the  E/H  field  sensors,  can  produce  a wave  im- 
pedance measurement  error  in  of  2.5:1.  (Effects 
of  a Grounded  Parasitic  Line  on  Earth  Conductivity 
isurements , ' I , . 14 , dated 

1966. ) 

: . E-i  Lnal  i ■ • < ( f roi  near- 

by telephone  cable  parasitic  element  when  the  common 
mode  impedance  is  less  than  200  ohms. 

c.  For  the  antenna  coni Lgui  ition  Pi  ire  14, 

it  would  seem  that  the  major  contamination  is  in 
the  E-field.  A long  wire  grounded  parasitic  near 
the  AMT  long  wire  E- field  sensor  would  have  the 
effect  of  moving  a ground  plane  closer  to  the  sens- 
ing antenna  (shielding),  thereby  reducing  the 
fields  measured  by  the  receiver.  The  H-field  sensor, 
on  the  other  hand,  is  not  as  sensitive  to  single- 
ended  long  wire  parasitics,  since  it  sees  only  a 
small  volume  of  the  area  over  which  the  parasitic 
H-field  exists.  For  the  case  of  a single  Long  wir< 
parasitic  line,  grounded  remotely  from  the  11-field 
sensor  locations,  the  H-field  contamination  would 
probably  be  less  than  10  percent. 

The  GTE  Schlumberger  (1972)  expansions  for  the  E/W  and 
N/S  antenna  lines  show  a relatively  consistent  anisotropic  . 
profile  in  both  directions  suggesting,  once  again,  an 
erroneous,  measured  5:1  a I between 

legs . 

...  Layeri  j Earth  Interpi  tation 

The  results  of  this  survey  are  in  general  agreement  with 
previous  GTE  Schlumberger  (1972)  studies  at  this  site  which 
modelled  the  following  layered  earth  structure: 
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E/W  LINE  N/S  LINE 


670  — m 

40 

m 

670  -m 

40 

m 

17000 

240 

m 

26700 

240 

m 

350 

340 

m 

270 

340 

m 

6800 

9035 

(4170) TM  (8910 )TE 


The  E/W  antenna  depth  sounding  curves  and  lateral  pro- 
filing can  be  interpreted  as  describing  a layered  earth 
structure  consisting  of  a low  resistivity  layer  at  the  sur- 
face followed  by  a high  resistivity  middle  layer  and  another 
low  resistivity  : iy<  r at  lepth.  ■ • 

probably  corresponds  to  the  upper  three  layers  of  the  Schlum- 
berger  data  where  the  high  resistivity  layer  is  not  detect- 
able. The  recent  AMT  results  at  79  Hz  (the  maximum)  should 
reflect  the  properties  of  the  lower  layer  determined  in  the 
Schlumberger  work. 


N/S 

E/W 

AMT*  (maximum  value) 

8,910 

4 , 170 

Lower  Layers  Schlumberger 

9,035 

6 , 810 

*Two  polarizations  along  E/W  antenna  line. 

The  agreement  of  the  values  and  of  the  sense  of  the 
anisotropy  is  surprisingly  good.  This  data  and  the  full  an- 
tenna measurements  suggest  a shear  pattern  with  an  east-west 
strike  in  the  bedrock. 

The  depth  sounding  results  shown  in  Figure  15  represent 
the  average  resistivities  over  the  full  E/W  antenna  at  the 
sampled  intervals  for  two  polarizations. 

These  data  suggest  quite  clearly  the  presence  of  a four- 
layer  earth.  Because  the  top  layer  is  only  a few  tens  of 
meters  thick  and  quite  strongly  anisotropic  (if  these  data  are 
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taken  at  face  value),  a simple  three-layer  interpretation 
(a  hiqh  resistivity  layer  sandwiched  between  two  low  resis- 
tivity layers)  can  be  postulated.  Koziar  (Ph.D.  thesis,  m 
preparation)  shows  that  this  situation  can  be  inverted,  usinq 

the  peak  resis ilu<  md  thi  frequ<  tncy  it  thi 

The  followinq  conclusions  can  be  drawn: 


The  upper  layer  can  be  characterized  by  the  value 
of  h. /(  or  the  ratio  of  its  thickness  to  its 
resistivity 


h?  _ 

P 2 


.252 


( . max  f max; 


T/2 


> - k <-m 
f - kHz 


The  resistivity  of  the  middle  layer  cannot  be  de- 
termined precisely  but  must  be  max. 

The  thickness  of  the  middle  layer  is  qiven  by 

1 1/2 

h,  = 503 


p 

max 

f 

max 

d.  The  resistivity  of  the  bottom  layer  (.4)  . , and 

less  than  the  lowest  observed  value. 

This  leads  to  the  followinq  one-dimensional  interpreta- 
tion for  the  layered-earth  structure  at  the  WTF : 


Layer 

1 Surface  Stronqly  anisotropic  - a few  tens  of 

meters  thick 

2 p 1700  ohm-meter;  h2  600  meters 

3  ,3-6600  -m;  h 4 . 6 km 

4 p,  <-  1000  U-m 

The  full  antenna  AMT  measurements  alone  the  N/S  antenna 
add  a second  dimension  to  the  above  layered  earth  interpreta- 
tion . Th<  result!  t > • ti  • rtl  ti  I ex- 

hibit low  conductivities.  The  GTE  Sylvania  Schlumberqer  (1972) 
data,  combined  with  the  E/W  residual  plots  and  half-leq  measure- 
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merits,  suggest  a general  E/W  geologic  fracturing  in  the  bed- 
rock. A simplistic  model  would  consist  of  a shear  line  inter- 
secting a north  leg  1500  meters  from  the  center  and  extending 
to  the  west  leq  5600  meters  from  the  center.  Everything  north 
: 1 sheai  Lii  would  1 ve.  tro  eai 

the  center  of  the  array  could  be  as  high  as  2. 

5.2.4  ■ W : londuct 

The  results  from  the  October  AMT  lateral  profiling  and 
the  January  half  and  full  antenna  measurements  show  once 
again  thi  WTF  area  to  be  geoelectrically  complex  and  highly 
an i sot  roi  ic. 

The  survey  emphasis  was  on  conductivity  measurements  alonq 
the  WTF  E/W  antenna  line.  The  program  included  lateral  pro- 
filing from  20  measuring  stations,  half  and  full  antenna 
line  measurements.  The  analysis  showed  the  effective  con- 
ductivity for  the  E/W  antenna  line  at  80  Hx  to  be  2.6  x 10”  ’ 
mhos/meter.  The  results  of  this  survey  and  various  other 
WTF  studies  are  summarized  m Table  XI. 

Some  observations  can  be  made  from  the  results  of  this 
survey  and  past  WTF  studies.  These  are  as  f o 1 

• 

full  antenna  conductivities  are  in  jood  agreement, 
within  +8  percent.  Statistical  analysis  has  shown 
that  the  GTE  Sylvania  measured  E W antenna  con- 
ductivity data  has  the  lowest  standard  deviation 
of  all  other  independent  half  or  full  antenna 
measurements . 

• 

tive  and  numerical  averaging  of  many  short  electrode 
AMT  measurements  consistently  produce  values  lower 
than  long  line  measurements.  One  possible  explana- 
tion is  that  200’  AMT  electric  dipoles  do  very 
little  volume  a raging,  and  consequently  the  measure- 
ments along  the  E/W  vitenna  might  not  impac  the 
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high  conductive  regions  to  the  north  and  west. 

Eleven  kilometer  electrodes  would  certainly  see 
this  geologic  fracturing  much  sooner. 

• Sy  1 v.in 

previous  studies,  while  N and  n/s  values  seem  to 
be  in  error,  suggesting  parasitic  contamination . 

The  GTE  Sylvania  Schlumberger  (1972)  expansions  for 
the  E/W  and  N/S  antenna  lines  show  a relatively  con- 
sistent anisotropic  , profile  in  both  directions 
suggesting,  once  again,  an  erroneous,  measured  5:1 
anisotropy  contrast  between  the  N/S  legs. 
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APPENDIX  A 


MICHIGAN  DATA 


CONTENTS 


SECTION 


PAGE 


Log  Resistivity  vs.  Log  Frequency 
with  Polynomial  Fits 

North-South  Orientation  A-3 

East-West  Orientation  A-19 

Averaged  Data  A-35 

Psendosections 

North-South  Orientation A-52 

East-West  Orientation A-57 

Averaged  Data .A-62 

Residuals 

North-South  Orientation A-68 

East-West  Orientation A-78 

Averaged  Data A-88 

Anisotropy  (TE/TM  Radio) A-98 


Michigan  Station  Location A-110 


SECTION  A.  1 


MICHIGAN 


LOG  RESISTIVITY  VS.  LOG  FREQUENCY  WITH  POLYNOMIAL  FITS 
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SECTION  A. 2 


MICHIGAN 


PSEUDOSECTIONS 


• NORTH-SOUTH  ORIENTATION 

• EAST-WEST  ORIENTATION 
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SECTION  A.  3 

MICHIGAN 

RESIDUALS 

1 • NORTH-SOUTH  ORIENTATION 

• EAST-WEST  ORIENTATION 

• AVERAGED  DATA 
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SECTION  A. 4 


MICHIGAN 


ANISOTROPY 
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SECTION  A. 5 


MICHIGAN 


STATION  LOCATIONS 


SECTION  A. 5 


MICHIGAN  REGION  STATION  LOCATIONS 


LATITUDES  LONGITUDES* 


STATION  DESIGNATION 

DEG 

MINS 

DEG 

MINS 

Line  A 

S30 

46 

49 . 2 

271 

45.3 

S28 

46 

50.5 

271 

50.2 

S29 

46 

49 . 2 

271 

54 . 2 

M48 

46 

50 . 2 

271 

57.1 

S10 

46 

54.  1 

271 

57.8 

S11R 

46 

47.7 

272 

21.5 

Line  B 

S27 

46 

44 . 4 

271 

29.7 

M42 

46 

44 . 6 

271 

37.1 

S31 

46 

44 . 7 

272 

4 . 1 

S32 

46 

45.6 

272 

10.9 

S33 

46 

44 . 2 

272 

11.5 

S213 

46 

45.6 

272 

11.7 

S210 

46 

43.8 

272 

16 . 8 

Line  C 

S46 

46 

38.9 

271 

24 . 2 

S26 

46 

39.1 

271 

24 . 5 

S212 

46 

39.7 

271 

53.5 

S20  3 

46 

37.5 

271 

57.7 

N12 

46 

35.2 

2,2 

9 . 5 

S3  3A 

46 

37.3 

2i2 

10.3 

S3  4 

46 

37.3 

272 

12.8 

S211 

46 

36.8 

272 

21.8 

Line  D 

S202 

46 

33.1 

271 

37.9 

S43 

46 

32.5 

271 

44 . 3 

S204 

46 

32.5 

271 

59.8 

M35 

46 

34.9 

272 

0.1 

M94R 

46 

33.3 

272 

9.0 

S24 

46 

34.2 

272 

15.3 

Line  E 

S44 

46 

28.2 

271 

24.3 

S48 

46 

28.5 

271 

29.3 

N9 

46 

25.7 

271 

29.4 

For  true  longitude,  subtract  from  360°. 


LATITUDES 

LONGITUDES* 

STATION  DESIGNATION 

DEG 

MINS 

DEG 

MINS 

S42 

46 

28.7 

271 

38.2 

S45 

46 

25.0 

271 

45.0 

M89 

46 

29.7 

271 

51.4 

S205 

45 

27.5 

271 

54.9 

S206 

46 

27.7 

272 

5.6 

S207 

46 

27.2 

272 

17.6 

S208 

46 

28.2 

272 

27.7 

S3 

46 

26.7 

272 

33.5 

SI 

46 

27.0 

272 

35.0 

S4 

46 

25.7 

272 

37.7 

S4B 

46 

25.3 

272 

39.2 

S2 

46 

27.0 

272 

46.0 

Line  F 

S25 

46 

21.0 

271 

24.6 

S37 

46 

21.4 

271 

32.6 

S41 

46 

21.2 

271 

41 . 2 

S40R 

4 6 

21.3 

271 

45.7 

M105 

46 

22.7 

271 

54 . 5 

S13B 

46 

21.0 

271 

55.7 

M4 

46 

20.7 

272 

12.6 

S17 

46 

20 . 9 

272 

23 . 3 

S13AR 

46 

24.  1 

272 

26.2 

S13 

46 

22.1 

272 

29.9 

S5B 

46 

21.2 

272 

34.9 

S5 

46 

22.8 

272 

35.5 

S6R 

46 

19.3 

272 

38.4 

S7 

46 

20.2 

272 

43.0 

S12 

46 

21.2 

272 

44 . 3 

Line  G 

S47 

46 

16.1 

271 

25.8 

S3  8 

46 

16.2 

271 

34.7 

S39 

46 

16.2 

271 

39.4 

S87 

46 

16 . 1 

271 

46.1 

N3 

46 

16.8 

271 

51.9 

S49 

46 

17.2 

271 

54.6 

N16 

46 

15.3 

272 

2.9 

S21 

46 

16.2 

272 

8.4 

S16 

46 

16.7 

272 

15.3 

S50 

46 

15.8 

272 

17.4 

SI  5 

46 

15.5 

272 

26.3 

* For  true  longitude, 

subtract 

from  360°. 

A- 


LATITUDES 

LONGITUDES* 

STATION  DESIGNATION 

DEG 

MINS 

DEG 

MINS 

SI  4 

46 

15.6 

272 

32.2 

S8R 

46 

16.8 

272 

42.  3 

S9 

46 

16.7 

272 

48.3 

Line  H 

S71 

46 

11.2 

271 

23.8 

S73 

46 

10.9 

271 

29 . 3 

S84 

46 

10.8 

271 

40.5 

S67 

46 

11.5 

271 

45.2 

S81 

46 

10.7 

271 

50.2 

S72 

46 

11.0 

271 

53.1 

S69 

46 

11.1 

272 

1.1 

S68 

4r 

10.0 

272 

11.8 

S64 

4b 

11.4 

272 

22.0 

S51 

46 

12.4 

272 

23.8 

S52 

46 

11.5 

272 

28.0 

S53 

46 

12.2 

272 

33 . 8 

S54 

46 

10.3 

272 

40.4 

S55 

46 

8.5 

272 

43.7 

Line  I 

S86 

46 

7.1 

271 

31.9 

S75 

46 

5.9 

271 

46.1 

S74 

46 

7.1 

271 

51.1 

S79 

46 

5.7 

272 

0.5 

S77 

46 

5.9 

272 

6.2 

S70 

46 

6.9 

272 

9.6 

S61 

46 

6.7 

272 

14 . 3 

S62 

46 

6 . 7 

272 

20.2 

S60 

46 

6.5 

272 

26.9 

S59 

46 

6.1 

272 

30.7 

S57 

46 

5.8 

272 

39 . 8 

S56 

46 

5.7 

272 

44.1 

Line  J 

S85 

46 

1 . 4 

271 

40.7 

S82 

46 

1.8 

271 

47.0 

S83 

46 

1.7 

271 

52.0 

S80 

46 

0.9 

271 

58.8 

S130 

46 

1.8 

272 

6.3 

S7  8 

46 

1.4 

272 

8.8 
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For  true  longitude,  subtract  from  360°. 


LATITUDES 

LONGITUDES* 

STATION  DESIGNATION 

DEG 

MINS 

deg’ 

MINS 

S63 

46 

1.3 

272 

18.5 

S200 

46 

1.6 

272 

36.9 

S58 

46 

1.5 

272 

40.8 

* "or  true  longitude,  subtract  from  360°. 

■ 

Stations  not  used  (either 

repeats 

or  bad 

station 

measure 

ments)  : 

S 3 8 , S22 , Sll,  S11R,  S19,  S32R,  S21A, 
S4BR , S13A,  S6 , S8 , S40,  M94,  N12R 
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NUSC  SITES 

NUSC- 

SELECTED 

SITE 

LOCATIONS 

AVERAGE  LOG  p 
(75  Hz) 

LATITUDES 
DEG  MINS 

LONGITUDES 
DEG  MINS 

N2 

46 

36.9 

87 

36.1 

3.25 

N3 

46 

18.1 

88 

8.8 

2.95 

N4 

46 

11.5 

87 

58.9 

3.10 

N5 

46 

23.2 

87 

8.4 

1 . 77 

N6 

46 

14.3 

87 

15.7 

3.68 

N7 

46 

12.1 

87 

28.8 

3.64 

N8 

46 

13.5 

88 

21.1 

3.11 

N9 

46 

25.7 

88 

30.5 

2.54 

N10 

46 

39.0 

88 

17.8 

3.28 

Nil 

46 

54 . 2 

88 

5.0 

2.46 

N12 

46 

34.3 

87 

50.4 

4 . 44 

N14 

46 

20.3 

88 

19.1 

3.55 

N15 

46 

30.5 

88 

16.7 

3.77 

N16 

46 

14.8 

87 

54.4 

3.70 

N28 

46 

23.2 

88 

39.2 

2.39 

N32 

46 

2.4  ' 

87 

41.9 

2.36 

N33 

46 

7 . 1 

87 

48.5 

3.77 
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